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Abstract We study packings of n hard spheres of equal radius in the d-dimensional 
unit cube. We present a nonsmooth function whose local extrema are the radii 
of jammed packings (where no subset of spheres can be moved keeping all others 
fixed) and show that for a fixed number of spheres there are only finitely many radii 
of such jammed configurations. We propose an algorithm for the maximization of 
this maximal radius function and present examples for 5-8 disks in the unit 
square and 4-6 spheres in the unit cube. The method allows straightforward 
generalization to packings of spheres in other compact containers. 
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1 Introduction 

Originally introduced as models for monatomic liquids, granular media and glasses, 
random dense sphere packings have been of high interest to experimental and the- 
oretical physicists as well as mathematicians for more than 50 years. The classical 
experiments of Bernal [2,3] and Scott [28] from 1960 were to pour small balls of 
equal radius into a container and shake the container for the packing to densify. 
Although this is a parameter-dependent protocol (the speed of pouring and the 
shaking amplitude and frequency), it was observed that there exists a "universal", 
highly reproducible packing fraction of « 0.63 — 0.64 for disordered monodisperse 
spheres [17]. Later, Lubachevsky and Stillinger [22] introduced the "stochastic bil- 
liard" , a computational tool to create dense sphere packings. There, a fixed number 
of spheres move within a container while their radii increase at a rate common to 
all. The spheres collide elastically with each other and with the walls of the con- 
tainer if such are present. At every collision, the kinetic energy of the collision 
partners is increased to ensure their separation. This procedure is stopped once 
either the collision frequency or the pressure (the sum of the squared velocities) 
exceed a certain threshold. This is also a parameter dependent protocol, in that 
the radius growth rate has to be set by the user. Very slow growth rates result in 
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higher packing fractions. However, the spheres are then arranged in local patches 
of crystalline, and thus regular substructures. Applications of random dense sphere 
packings are manifold, from modeling the burning of solid rocket propellant [18], 
predicting the drug release kinetics of matrix tablets in pharmaceutical science [1] 
to phase transitions in statistical mechanics [10,20], to name only a few. 

Torquato et al. [29] argued that the concept of random close packing of spheres 
is ill defined. Clearly, the properties of density and randomness are at variance 
with each other and an increase in one must come at the cost of a decrease of the 
other. The authors of [29] propose to balance an increase in the packing fraction (f> 
with an increase of a suitably chosen "order parameter" ip. They suggest that there 
are "maximally random jammed" (MRJ) states, namely jammed sphere packings 
that minimize that order parameter. Following Torquato et al. [29,12], a packing 
of finitely many spheres in a container is (collectively) jammed if no subfamily of 
spheres can be displaced continuously while fixing the positions of all other spheres. 
The present paper is motivated by the question for which radii there are jammed 
configurations of n spheres in the d-dimensional unit cube. Many authors have 
studied densest packings of 50 - 100 spheres in squares, triangles, disks and the 
cube, for certain numbers of spheres also with computer-aided proofs of optimality, 
see for example [4,9,13,16,24,25,26,27] and the references therein. Nurmela and 
Ostergaard [25] achieved this by minimizing the (smooth) discrete Riesz s-energy 
of a configuration x of sphere centers x^ (for more details on the notation see 
definition (1) below) 



^(x) = ^|x i -x,-| s 

i<j 

in the limit s — >■ oo, when this energy approaches the nearest neighbor distances 
[14], see also the recent work by Carlsson et al. [6]. Rather than seeking only 
densest packings, here we emphasize finding as many jammed packings as possible 
and work directly with a nonsmooth function to be optimized. 

In this paper we take up the idea from [27] of the maximal radius function G 
(see Equation (2) below) for finite configurations and show that jammed packings 
correspond to its local maxima. We prove that for every finite number of spheres 
in [0, l] d there are only finitely many radii at which jamming can occur. Then we 
apply techniques from nonsmooth analysis [7, 23] to optimize G. Our algorithm 
is based on a line search and has similarities to the gradient sampling algorithm 
by Burke et al. [5]. To find the direction of the line search we exploit the special 
structure of the objective function whose gradients are readily calculated. This op- 
timization procedure is repeated many times with random starting points. With 
this parameter-free approach we are able to find many jammed packings for up to 
eight disks in the square and six spheres in the cube together with the minimal 
polynomials for the corresponding radii (thereby proving their existence). To the 
best of our knowledge, many of these configurations have not been reported else- 
where in the literature. Our method can easily be generalized to sphere packings 
in arbitrary convex polytopes, within a sphere or to spherical caps on a sphere. 
The paper ends with a discussion of open problems for future research. 
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2 Jammed sphere packings and the maximal radius function 

Throughout the paper, | • j denotes the Eucliden norm. Let the number of spheres 
n > 2 be fixed and let 

x = (xi, ... ,x„) = (xu, . . . ,x ld , . . . ,x nl , . . . ,x nd ) (1) 

be a collection of n points Xj £ M. d , d > 2. We say that x is r- admissible for r > 
if |xj — Xj > 2r for all i, j = 1, . . . ,n with i ^ j and r < x ik < 1 — r for all 
i = 1, . . . , n and k = 1, . . . , d. The latter conditions can be easily generalized to 
convex polytopes bounded by half-spaces with known normal vectors. Notice that 
for any such point x, a permutation of the entries Xj gives another point x' that 
determines the same geometric configuration of spheres. In an even stronger sense, 
one can identify configurations that are obtained from each other by an application 
of a symmetry of the cube. 

Denote by M d (r) C M. dn the set of all r-admissible configurations of n dis- 
tinct points in [0, l] d . When r = 0, then M„(0) is path-connected, since any 
two components Xj and Xj can exchange their places along continuous curves, 
each avoiding all other particles and the faces of the cube. If r\ < r<2, then 
Mn{r2) C Mn(ri). Moreover, for every d there is a sequence v n with lim r n — 

71— >00 

such that M d (r n ) = (large numbers of spheres imply small radii). As r increases, 
path-connected components split off (i.e. any two points in the same component 
can be joined by a continuous curve) at certain times and disappear at others 
and so the topology of M d (r) changes. We say that r is critical, if there exists a 
£ > such that for the numbers of connected components (denoted by /3o) we 
have fio(Mn(s)) j= Po{M%(t)) for every r-E<s<r<t<r + e. In particular, if 
Po{Mn(s)) > Pa{M d (t)), we speak of a disappearance and we call the set 

C(r) = M d (r) \ p| Mi(t) 

t>r 

a critical set (this is the union of all disappearing connected components at the 
critical radius r). 

For r > 0, M d {r) is clearly closed. If x € int M d {r) and r > 0, then no two 
spheres centered at the points in x touch and there exists a e > such that 
x G intM d (r') for every r' S [r,r + e). Thus C(r) at a critical radius has empty 
interior. We call a point x G C(r) a partially jammed configuration of spheres and 
an isolated point of C(r) a fully jammed configuration of spheres. The fact that 
intC(r) = excludes that the spheres can be displaced continuously so that 
eventually they all loose contact with the walls and each other and the radius can 
be increased. However, it is possible that a partially jammed configuration leaves 
room for an unconstrained sphere (a "rattler", see Figures 3 and 4) 1 . Clearly, there 
is always a largest radius r for which M d (r) is not empty. 

Jammed configurations can be characterized as local maxima of the maximal 
radius function that we now introduce. With the notation as in (1), let 

<Pij(x) = 2 , i = l,...,n- 1, j = i + ...,n, 

■0ifeW=a!,- fc) 4< lk {x) = 1 - x lk , i = l,...,n, k = 



1 Figures and Tables are collected at the end of this manuscript. 
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and denote the set of these N(n,d) := n 2 + 2nd functions by T. Define 
G : [0, l] nd — > [0,oo) by the lower envelope 

G(x) = mm{/(x)}. (2) 

This is the maximal r such that a sphere of radius r can be centered at every entry 
of the n-tuple x = (x.%, ... , x n ) without violating any of the other spheres or the 
walls of the unit cube. Throughout, we fix an index set £ = {1,2,..., N(n, d)} and 
relabel the functions so that T = {/; : / G £}. The gradients of the component 
functions are 



V0,j(x) = — — 1 I . . . ,0, . . . ,x ik - Xjk, . . . ,0, . . . ,xjk - x ik , . . . , 0, 

Z|Xj Xj 



where the first stretch begins at the entry (i — l)d and the second begins at the 
entry (j — l)d. The gradients of the (j>ij exist wherever x^ ^ x^ and are non-zero 
there. Trivially, 

V^ife(x) = 5((i_i)rf + fe,i), Vip lk {x) = 

where S k i denotes the Kronecker symbol (1 if k = I and otherwise). 

Lemma 1 7/x* is a local maximum of G, then x* is a partially jammed configuration. 
If x* is a strict local maximum of G, then it is fully jammed. 

Proof. We have that M^(r) = G _1 ([r, 00)). If there exists a 8 > such that 
G(x) < G(x*) for all x x* with |x — x*| < 8, then none of these x lies in 
M^(G(x*) + e) for any e > 0. □ 
We use methods from nonsmooth optimization to find the maxima of G, see [7, 
23] for references. For a A"-Lipschitz continuous function, the generalized directional 
derivative at x in the direction v is 

f o, v ,. /(y + tv) - /(y) 
/ (x;v)= hmsup — f 

y->x,t\,0 t 

this is bounded from above by A'|v|. The generalized gradient or the Clarke subdif- 
ferential is the nonempty set 

a/(x) = {( £ R m : /°(x; v) > C ■ v for all v G R m } . 

If / has an extremum at x, then G <9/(x), [23, Theorem 3.2.5]. For a function G 
defined by a minimum selection of functions indexed by a set C as in Equation (2) 
let 

£(x)={(e£:/,(x)=G( x )} (3) 

be the indices of the functions that realize the common minimum; these are called 
the active indices. Then the subdifferential of G at a point x is 



<9G(x) = conv{V//(x) : / G £(x)}, 
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the (closed) convex hull of the active gradients, this follows for example from [7, 
Theorem 2.5.1]. Moreover, the generalized directional derivative is 

G°(x;v) = min{V//(x) • v : I 6 £(x)}. 



Proposition 1 For every number of spheres n and dimension d, the function G has 
finitely many local maxima. 

Proof. This follows from the fact that every subset of J" can be active at one local 
maximum at most. Assume to the contrary that there are local maxima x and 
x' with £(x) < £(x') and G(x) < G(x'). As all functions participating in T are 
convex, the sublevel sets /j _1 ([0, G(x)]) are convex. Since the higher point x' is a 
local maximum and hence G conv{V//(x') : I £ £(x )}, there must be at least 
one gradient that satisfies V/;* (x') • (x' — x) < 0. Because of the convexity of the 
sublevel sets of /;» this would imply G(x) > G(x'), a contradiction. □ 
In the following we describe a iterative procedure to find local maxima of the 
maximal radius function G. An initial choice x is selected for example randomly 
from a uniform distribution in [0,1] . Given the iteration x fe , let J(x fc ) be the 
matrix whose rows are the active gradients V/;(x fe ) for all I € £(x fc ). We seek a 
direction £ fe that solves the minimization problem 

i k e argmin|J(x' c )£-l| 2 , (4) 
S 

where 1 is a vector with |£(x fc )| entries 1. This implies that the active functions 
all increase infinitesimally at the same rate. With such a £ in hand, we begin a 
line search in that direction. A triple < t\ < t 2 < t% is a bracket of a directional 
maximum, if 

G(x fc + ti£ fc ) < G(x fc + t 2 £ k ) > G(x ft + t 3 £ k ), (5) 

and at least one of the inequalities is strict. If the width of the bracket t 3 — t\ is 
sufficiently small, then we set x fe+1 = x fc + t 2 £ k ■ 

There are three issues that need to be addressed in the numerical implementa- 
tion of this algorithm. First, due to rounding errors we need to relax the equality 
requirement in the definition (3) of the set of active functions C. This is done by 
replacing (3) by 

£e(x) = {Ze£ : /;(x) <G(x)+e}, 

where the tolerance e can be adjusted if need be. For example, if active functions 
are lost after an iteration, then it is helpful to increase e. The second is a detection 
and treatment of saddle points. At the beginning of the line search procedure we 
set t\ = and t 2 to the machine precision. If instead of the first inequality in 
(5) we have the reverse strict inequality, then we perturb x fc by a random vector 
whose components are normally distributed with standard deviation 10 -3 and 
recalculate the direction of increase £ fc accordingly. Finally, the iteration terminates 
if the optimal function value from the minimization problem (4) exceeds a certain 
threshold to be set by the user or a maximum number of iterations has been 
reached. Once an approximate local maximum has been found, it can be refined 
by solving a system of linear-quadratic equations determined by the expected 
contact graph (using e.g. Mathematica). See Figure 1 for an example. 
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3 Examples 

The algorithm described in Section 2 has been implemented in the open source 
package scilab [11] and together with a Mathematica notebook will be made 
available on the author's website. 

The minimal example for multiple jammed configurations of different radii is 
that of five disks in a square see Figure 2 (top row). In the first two cases, the 
critical set consists of 4 • 5! = 480 isolated points (one of the four corners is dis- 
tinguished). In the last case, the critical set consists of 4! = 24 isolated points 
only, since this arrangement has a greater degree of symmetry. We can identify 
the cascade of radii at which the number of connected components of Mf(r) in- 
creases (Figure 2, bottom row). We know various jammed configurations of six to 
eight disks (see Figures 3, 4 and 5). We also find the minimal polynomials whose 
smallest positive roots are the maximal radii with the help of the Mathematica 
function GroebnerBasis, see Tables 1 and 2. In cases where it is not obvious, a con- 
figuration can be tested for full jamming by Connelly's criterion [8, Equation (1)]. 
Namely, if there is an infinitesimal motion x' = (x 1; ... , x^) of the configuration 
x = (xi, . . . ,x„) that satisfies 

(xi - Xj) • (x- - x'j) > 

for every pair of touching disks, then the motion should vanish, i.e. x' = 0, 
for the configuration to be fully jammed. 

We repeat the maximization procedure described in Section 2 about 10 4 times 
with initial choices x° distributed uniformly in [0, l] nd . Locally maximal config- 
urations are compared to a list of known maxima (listed with all their images 
under symmetries of the square) by successive minimization of distances between 
the points of the new candidate and the known configuration. Although the maxi- 
mization procedure is not entirely deterministic due to the randomized symmetry 
breaking, this gives a rough measure of the relative sizes of the basins of attraction 
of the different maxima. It is a general observation that denser configurations are 
reached more frequently, see the right panel of Figure 6. 

Contact graphs of packings of four to six spheres in [0, 1] are given in Figures 
7 - 9 and Tables 3 - 5. 

We end this section with some test configurations that allow other researchers 
to benchmark their methods. It has to be stated again that our algorithm has a 
random component and the path of the algorithm to the local maximum cannot 
be reproduced exactly. The first is a configuration of seven disks in the square, the 
second is a configuration of five spheres in the cube, and the third is a configuration 
of five disks in an equilateral triangle of side length 2, 

x = (0.8, 0.1; 0.9, 0.8; 0.1, 0.9; 0.9, 0.2; 0.6, 0.3; 0.1, 0.5; 0.2, 0.2), 
x = (0.1,0.9,0.8; 0.2,0.1,0.9; 0.9,0.2,0.3; 0.6,0.1,0.4; 0.5,0.3,0.2), 
x = (0.6, 0.3; 1.0, 0.4; 1.5, 0.2; 0.9, 1.2; 1.2, 1.3). 

Figures 10 - 12 show the results of runs with these initial configurations. 
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4 Discussion 

The above methodology may be applied to other geometrical setups, for example 
spheres in a regular simplex or spherical caps on the unit sphere S d (also known as 
Tammes problem). It is also possible to have multiple radius classes with fixed radius 
ratios. As stated in the introduction, the goal of the work is not to find the global 
maximum of the maximal radius function G, but rather as many local maxima 
as possible. Unfortunately, there is no way of knowing how many local maximal 
values there are, nor to know the "basins of attraction" of the local maxima for the 
given (or any other) search algorithm. There is a choice of alternative algorithms 
for a single optimization of the function G, for example the gradient sampling 
algorithm [5] or derivative-free methods [15]. One also has to study alternative 
choices for the distribution of initial points (random or not), in order to find as 
many local maxima as possible. This touches upon the "curse of dimensionality" 
which is a well-known problem in numerical integration (for example over [0, l] d ), 
see [19] for a recent review. 

It is possible to improve the performance of the search algorithm by exploiting 
the special structure of the objective function G further. Among the 0(n 2 ) pair- 
wise distances between sphere centers, many need not to be computed for finding 
the minimum. For example, one can subdivide the unit cube into n cubes of equal 
side lengths and only compute those distances between sphere centers in adja- 
cent subregions [22]. This introduces some overhead, since membership lists of the 
subregions must be maintained. 

Although the data in Figures 2, 3, 4, 5 and 8 are sparse, they allow to make 
some interesting observations. First we observe that no contact graph appears at 
more than one critical value. This leads to interesting combinatorial problems. Is 
it possible to give upper and lower bounds on the number of local maximal values 
of a minimum selection of convex functions as in Equation (2)? A very coarse 
upper bound would be the number of subsets of the index set C with nd + 1 and 
more elements. Since the functions xpif. and ip lk are complementary to each other, 
at most one of these can be active and so a slightly reduced upper bound for 
minimally determined maxima would be 

/ n(n - l)/2 + nd\ 
\ nd+1 ) ' 

A starting point of this work was the question which packing fractions are 
realizable by strictly jammed packings. For a jammed packing x* of n spheres in 
[0, l] d , the packing fraction is 

<t>n,d = noj d G(x.*) d , 

where uj d is the volume of the d-dimensional unit ball. Here only few general things 
can be said, apart from the numerical area spectra in the right panel of Figure 6. 
In every dimension, there is at least one packing fraction that is realized infinitely 
many times, namely which can be achieved by any number of n d spheres. 
Other repeated packing patterns are known for disks [21,25]. 

To end this paper with an open question, is there an open interval in [0, 1] of 
packing fractions that are realizable as limits of finite packings as n — > oo? 
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Figures and Tables 




Fig. 1 From the approximately locally maximal configuration, it is possible to construct the 
expected contact graph by adding the dashed line. This leads to the system of equations on 
the right. Solving this system with Mathematica gives r = tt,/^,~m ~ 0.254333095 as the 

2[ 3+2 v 2 ) 

maximal value. 
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Tf (11) T| (11) 1* (12) 




Si (10) Si (11) Si (10) 




Fig. 2 fTop row) Five fully jammed disks in a square (with increasing radii from left to right). 
The straight lines indicate the contact graphs. The configuration Tj 5 allows up to three more 
disks of the same size to be added. Numbers in parentheses indicate the number of functions 
that realize the common radius. (Bottom row) Critical points at which connected components 
of M| (r) split off. The numerical ordering of the critical values is Sf < Tj 3 < 5| < 5| < T| < 




Fig. 3 Partially and fully jammed configurations of six disks in a square with increasing radii 
from left to right and top to bottom. Configuration T| allows another disk of the same size to 
be added. 
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Fig. 4 Partially and fully jammed configurations of seven disks and splitting events with 
increasing radii from left to right and top to bottom. The thin solid lines in configuration Tg 
indicate that it is indeed jammed. Configuration Tg also exists with the two disks as rattlers. 




Fig. 5 Partially and fully jammed configurations 
to right and top to bottom. 



of eight disks with increasing radii from left 
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Fig. 7 Contact graphs of packings of four spheres in [0, l] 3 with increasing radii from left to 
right and top to bottom. Note that the two configurations labeled Tj have identical contact 
graphs, but different arrangements of the spheres in the cube. 
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Fig. 9 Contact graphs of packings of six spheres in [0, l] 3 with increasing radii from left to 
right and top to bottom. 
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Fig. 10 (Left) An initial configuration of seven disks in the unit square and a local maximum 
provided by the search algorithm. (Right) The corresponding value of the function G over 200 
iterations of the line search procedure. 




Fig. 11 (Left) An initial configuration of five spheres in the unit cube and a local maximum 
provided by the search algorithm. (Right) The corresponding value of the function G over 200 
iterations of the line search procedure. 




Fig. 12 (Left) An initial configuration of five disks in the an equilateral triangle of side lenght 
2 and a local maximum provided by the search algorithm. (Right) The corresponding value of 
the function G over 200 iterations of the line search procedure. 
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2 - 38* + 293* 2 - 1180i 3 + 2632t 4 
1 - 8* + 14i 2 
-13 + 52t + 92t 2 



3104t J 



■ 1524t° 



37 - 472t + 1872t 2 - 2368t 3 + 832t 4 

4 - 64t + 368t 2 - 896t 3 + 769* 4 
1 - 16* + 92t 2 - 224t 3 + 193t 4 
13- 112t + 208t 2 

25 - 440* + 4656t 2 - 34176i 3 + 144896t 4 - 311296* 5 + 262144* 6 
-1 + 8* - 8t 2 - 32t 3 + 32t 4 

5 - 138* + 1650t 2 - 11280* 3 + 49316* 4 - 147424* 5 + 316132* 6 
-495168* 7 + 547136t 8 - 400512* 9 + 196096t 10 - 58368* 11 + 9216t 12 
481 - 5680t + 22856t 2 - 40416* 3 + 32016* 1 - 9344t 5 + 256t 6 

1 - 8t + 13t 2 



see Table 2 

9 - 342i + 5850* 2 - 59538i 3 + 400966* 4 - 1876622* 5 + 6223518* 6 
-14613362* 7 + 23977405* 8 - 27600528* 9 + 25173396* 10 
-22602672* 11 + 14091408* 12 

20 - 300* + 1621t 2 - 41524 3 + 5832t 4 - 4320t 5 + 1296i 6 
1 - 50* + 1134t 2 - 15474* 3 + 142078* 4 - 929874t 5 + 4481482t 6 
-16190694* 7 + 44168181* 8 - 90780660* 9 + 138781596* 10 
-153497992* 11 + 116495620* 12 - 54495056* 13 + 11899792* 14 
see Table 2 

1 - 24* + 248t 2 - 1472t 3 + 5640t 4 - 146244 s + 25248t 6 

-26240* 7 + 12304t 8 

37 - 472t + 1872t 2 - 2368t 3 + 832* 4 



4 - 64* + 368t 2 - 896t 3 + 769t 4 

2 - 30t + 155t 2 - 310t 3 + 169t 4 

-73 + 584t - 344t 2 - 3296t 3 + 368t 4 

1 - 12* + 52t 2 - 120t 3 + 148* 4 

1 - 16* + 92t 2 - 224i 3 + 193i 4 

13- 112t + 208t 2 

-1 + 8* - 8i 2 - 32t 3 + 32t 4 



Table 1 Minimal polynomials of the locally maximal radii of disks in the unit square, which 
is always the first positive root. 
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A nonsmooth program for jamming hard spheres 
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value minimal polynomial 



1 - 16t -f 104?" 
5 - 20t + 4i 2 
5 - At + 2t 2 



352t 3 + 704i 4 - 1024t 5 + 1216i 6 - 768t 7 + 64t 8 



13 - 1041 + 240t 2 - 128t 3 + t 4 



1 - 8t + 20i 2 - 
1 - 16t + 104t 2 
9 - 36t + 4t 2 
5 - 20t + 4t 2 



16t 3 



Wt° 



352t 3 + 704t 4 - 1024t 5 + 1216i 6 - 768t 7 + 64t 8 



628849 - 35215544* + 92091652CK 2 - 14925909856* 3 
+167788673872* 4 - 1387008330496i 5 + 8720616629536* 6 
-42537931528576* 7 + 162698001135232* 8 - 489697674935296t 9 
+H56878795748352* 10 - 2126444318666752* 11 
+2994825062826240* 12 - 3163357034848256* 13 
+2454984222679040* 14 - 1420158379393024* 15 
+698664367800320* 16 - 324380948299776* 17 
+55492564615168* 18 + 84665382207488* 19 
-46125516718080i 20 - 14650005520384* 21 
+9630847598592* 22 + 3387487682560i 23 - 753326358528i 24 
-527777660928* 25 - 87879057408* 26 - 4697620480* 27 
+16777216* 28 

13 - lOAt + 240t 2 - 128t 3 + t 4 



1 - 8t + 20t 2 - 
1 - I6t + 96t 2 
see Table 4 
see Table 5 
9 - 36t + 4t 2 



16t 3 + t 4 

- 256t 3 + 544t 4 - 2304i 5 + 5632t 6 - 4096i 7 + 256i 8 



Table 3 Minimal polynomials of the locally maximal radii of four to six spheres in three 
dimensions. 
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Peter Hinow 



value minimal polynomial 

"Tf 17850625 - 3427320000* + 327817998528t 2 - 20821826923264* 3 
+987726417983232* 4 - 37313693603198976* 5 
+1168907639294074880t 6 - 31219963077836095488t 7 
+725429602195873475072* 8 - 14890754410484706631680t 9 
+273272074409891243776000t 10 - 4526854896736132756807680* 11 
+68222063007549667368159232t 12 
-941472366257163248946987008t 13 
+11962606417179071943721172992t 14 
-140607566392008233399332569088t 15 
+1534983670264170814477055475712* 16 
-15618180998282167274579315654656* 17 
+148566699694213641005963068178432t 18 
-1324815392234110934918500019863552t 19 
+11101518878460160647263724506316800t 20 
-87608344468411144802040068910874624t 21 
+652374475492830725841694709689876480* 22 
-4592111614294706899389248431215083520t 23 
+30605659141869032519917088312592695296t 24 
-193427325865007346716195666391281434624t 25 
+1160816883779208050449949294650825113600t 26 
-6623654558628870509077400504674421833728* 27 
+35977991812578058500230628908181005271040< 28 
-186234401198438304242849849705495082827776t 29 
+919617920841624229559217715864229182963712t 30 
-4335933763403022506485823748571063052664832* 31 
+19536660186433782946138183434741460248821760* 32 
-84185159953565880143925521370127290339426304t 33 
+347153909291483094403651698938794743133372416* 34 
-1370726446384907196824324909361586404962861056i 35 
+5184639756910668175364231224718241744740679680t 36 
-18792199494194327747086885128090261375386910720t 37 
+65288329365056041796420742462877146325874376704t 38 
-217448630736569356568147184944646006522920304640t 39 
+694326625113846229591979032844468100187769798656t 40 
-2125402526658739058303714747630549925702678872064i 41 
+6236455238395122619884513811694349682546780930048t 42 
-17537647049426855545851198344255050069316202397696t 43 
+47253288368880129538720061055878544735369533849600t 44 
-121951401628459296647005792037531770690788067901440i 45 
+301357940570851601462385042150371731014485373091840t 46 
-712770888641580130614333642358451074651177636331520t 47 
+1612891265948281610269936395233794346617846429122560t 48 
-3490216403945726551834895659369234775060744696233984t 49 
+7219138516017948554073970044709643998449131791056896i 50 
-14265527844268055655822641834843546109671242475765760t 51 
+26917438040439551303328998818949515655121390821965824t 52 
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-48471547386879781831770952534709393150828442619478016* 53 

+83253087740653914752431225193771295139744614707101696* 54 

-136306032151071672536241697958629229037463749324177408t 55 

+212597617702783995411019367614122565664623794996641792t 56 

-315676397473038371399275005328607494391181053884628992t 57 

+445923436728234832244994572550328731272800693297086464t 58 

-598806842242809930669704195963713773883830547081330688i 59 

+763786818077750680401448080069900372795670519535894528t 60 

-92457228842077274864420901913833164105113889897709568(K 61 

+1061175825386612849942695930583393999857612398744567808t 62 

-1153644758729836865896398564033754917894620519247904768i 63 

+1186634969364566686182010652820942165674057471044354048t 64 

-1153455507584389481769581120122977784941191436716474368i 65 

+1058161638153990547563625448826876248032168945639227392i 66 

-914835663828729166165471137790642289331085778821840896t 67 

+744187348818469308993553344010646859561311884061704192i 68 

-568596753514594051724592403539423821769311269813223424t 69 

+407248835742964482814890145907615762308072363154997248i 70 

-272836210124234733862768949209622314217163293068361728t 71 

+170557758832305386465277389754206931191924956380266496t 72 

-99215033862300201842972147687337951993209296674357248t 73 

+53539566761294723788199878504166278615696624851091456* 74 

-26707644285492858965320763857969512827073056556974080t 75 

+12266288629855639900270649057381011218333129081094144t 76 

-5162943193820983687414910787305081452370186372907008i 77 

+1980868601852602989513421247035536681150408206319616t 78 

-688426326927744918898952117680288113777061687984128i 79 

+215116218168063036670267319908135566372063166857216* 80 

-59901796639027703320498208900405319734826239000576t 81 

+14705581801215668409771970575321818506353256169472t 82 

-3141086480295116234155128029356961578965813690368t 83 

+574378710574108291290991339995528442291169525760* 84 

-88169677048730063219192101145477751204231512064t 85 

+11117813722503308733184750051934924875071750144t 86 

-1135155931044179239223315832239384127376195584t 87 

+97001378040253342688399202622900447232393216t 88 

-8255498437868897949003417002125119706890240* 89 

+870913736231297124613164904984002472968192* 90 

-98453434885687096630438107630143209996288i 91 

+9179886398371767408215264742671120334848* 92 

-651588097346297563395669662197706915840* 93 

+35581606297548954230865050104772952064t 94 

-1423200492877425383193585529990938624i 95 

+29623222488105672354653249885700096t 96 



Table 5 Minimal polynomial of configuration TP. 



